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ABSTRACT 


Starting  with  the  Hamiltonian  for  a  linear  harmonic  chain  of  2N 
particles  of  mass  m  and  one  of  mass  M  we  have  carried  out  numerical  calcula¬ 
tions  for  the  momentum  autocorrelation  function  of  the  mass  defect  particle 
for  chains  with  finite  number  N  of  mass  points  and  for  non-zero  values  of  the 
mass  ratio  p  s  ni/M.  These  results  have  been  compared  with  the  well  known  expo¬ 
nential  relaxation  of  the  momentum  autocorrelation  function  which  is  found  to  be  the 
rigorous  result  when  passing  to  the  thermodynamic  and  weak  coupling  limit.  In 
these  limits  the  dynamics  of  the  mass  defect  particle  is  exactly  described  by 
•  Fokkcr-Planck  equation,  i.c.  a  stochastic  equation  of  motion.  Wc  have  shown 
that  to  an  excellent  approximation  an  exponential  relaxation  of  the  momentum 
autocorrelation  function  is  obtained  for  mass  ratios  as  high  as  m  *  0.1  and 
for  chains  with  only  fifty  particles.  Thus,  for  the  harmonic  chain  considered 
here,  the  stochastic  aquations  of  motion  can  be  applied  to  a  very  good 
approximation  far  outside  the  usually  imposed  thermodynamic  and  weak  coupling 
limits. 


Kay  Words:  Stochastic  Equations,  Fokkor-I’lanck  Equation,  Thermodynamic  Limit, 
Weak  Coupling  Limit,  Momentum  Autocorrelation  Function,  Linear 
Harmonic  Chain. 
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1.  INTRODUCTION 


In  the  derivation  of  the  master  equation  of  non-equilibrium  statistical 
mechanics  from  the  Liouville  equation  two  limiting  processes  are  always  invoked. ^ 
First,  one  performs  calculations  in  the  thermodynamic  limit  in  which  the  size 
of  the  system  and  the  number  of  particles  in  the  system  arc  allowed  to  increase 
without  bound  in  such  a  way  that  the  concentration  (number  per  volume)  of  the 
particles  remains  finite.  The  other  limit  is  the  weak  coupling  limit  where 
some  interaction  parameter  is  allowed  to  approach  zero  while  the  time  approaches 
infinity  such  that  their  product  is  a  constant.  The  passage  to  these  limits 
is  essential  in  the  rigorous  derivation  of  the  master  equation  from  the 
Liouville  equation, ^ 

From  a  practical  point  of  view,  however,  all  physically  realizable 
systems  are  finite,  the  interaction  parameter  docs  not  go  to  zero  and  it  is  too 
time  consuming  to  make  observations  as  time  goes  to  infinity.  The  question 
thus  arises  as  to  just  what  error  is  involved  in  applying  stochastic  equations 
to  finite  systems  with  finite  strength  interactions.  It  is  to  this  question 
that  this  paper  is  addressed. 

Our  model  system  for  this  study  is  the  ubiquitous  linear  harmonic  chain 
and  the  specific  property  to  be  investigated  is  the  momentum  autocorrelation 
function  of  a  mass  defeat  particle  in  the  chain,  h'e  have  chosen  this  model 
since  it  is  possible  to  carry  out  exact  dynamic  calculations  for  the  auto* 
correlation  function  for  finite  chains  with  finite  interaction  forces.  These 
can  then  be  compared  with  known  analytic  stochastic  results  obtained  in  the 
thermodynamic  (infinite  chain)  and  weak  coupling  (infinite  mass  defect)  limits. ^ 
Since  the  ultimate  aim  of  any  dynamic  theory  is  the  calculation  of  some  average 
time  dependent  quantity,  the  examination  of  the  momentum  autocorrelation 
function  is  an  appropriate  test  for  the  validity  of  the  stochastic  equations. 
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(3) 

It  is  known'  1  that  the  dynamical  behavior  of  local  properties  in  a 
finite  chain  of  N  oscillators  with  short  range  forces  is  essentially  identical 
with  the  N-**>  results  over  a  range  of  time  which  is  proportional  to  the  size  N  of 
the  system.  We  present  a  numerical  study  of  this  situation  which  clearly 
demonstrates  this  point.  We  find  that  for  N  50  the  computer  results  on  the 

relaxation  of  the  momentum  autocorrelation  functions  are  indistinguishable 
from  the  thermodynamic  limit  results  over  a  range  of  time  in  which 
the  autocorrelation  function  has  relaxed  essentially  to  its  zero  value.  This 
certainly  indicates,  at  least  for  the  model  studied  here,  that  the  stochastic 
equations  arc  valid  for  finite  systems  far  removed  from  the  thermodynamic  limit. 

We  have  also  examined  numerically  the  weak  coupling  limit.  It  is 
known'  '  that  in  a  one  dimensional  harmonic  lattice  with  point  masses  m  and  an 
impurity  of  mass  M,  the  heavy  particle  undergoes  Brownian  motion  in  the  weak 
coupling  limit  w  m  m/M  **0,  t  •*  •,  pt«c.  This  implies  that  in  this  limit  the 
momentum  autocorrelation  function  decays  exponentially.  Our  nuamrical 
computer  study  shows  that  this  exponential  decay  of  the  momentum  autocorrelation 
function  is  also  obtained  to  within  a  very  close  approximation  for  non-zero 
values  of  u  and  for  finite  time*.  Thus,  excellent  agreement  with  exponential 
relaxation  is  already  found  for  w  ■  0.1,  i.e.  for  M  •  10m.  Clearly,  the 
stochastic  results  arc  again  la  excellent  agreement  with  the  exact  dynamic 
results  far  from  the  weak  coupling  limit. 

He  have  also  studied  the  infinite  chain  with  long  rango  interactions. 

We  have  shown  that  one  obtains  am  excellent  approximation  to  the  exponential 
relaxation  of  the  autocorrelation  function  of  the  heavy  impurity  M  for  w  >  0, 
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and  0  <  y  <  I#  where  y  is  an  interaction  parameter  such  that  y-»0  corresponds 
to  nearest  neighbor  interactions  and  y-»l  corresponds  to  interactions  of  zero 
strength  and  infinite  extent. 

In  Section  II  we  introduce  the  model  Hamiltonian  and  indicate  the 
various  representations  from  which  one  can  compute  the  momentum  autocorrela¬ 
tion  function  of  the  impurity.  In  Section  III  we  discuss  the  autocorrelation 
function  for  finite  systems  and  for  finite  interaction  strength  for  chains 
with  nearest  neighbor  interactions.  In  Section  IV  we  introduce  the  chain  with 
long  range  interactions  and  discuss  the  new  features  which  occur  in  this  system. 
In  Section  V  we  present  a  summary  and  discussion  of  our  results. 
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11.  HAMILTONIAN  AND  RliPRI.SLNTATlON  OF  AUTOCORRLLATI ON  FUNCTION 

The  Hamiltonian  appropriate  to  a  chain  of  2N  particles  of  mass  m  and 
one  of  mass  M  bound  by  harmonic  nearest  neighbor  forces  with  cyclic  boundary 
conditions  is 


a2 

p 


#A  A  .2 

<VVi> 


+ 


a 

2 


(2.1) 


Here  the  are  the  momenta  and  positions  of  the  2N  equal  mass  particle, 

P  and  Q  arc  the  momentum  and  position  of  the  impurity  particle  and  a  is  the 
force  constant.  If  one  partially  diagonalizes  this  Hamiltonian  according  to 
the  transformations 
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whore  w  is  the  fundamental  frequency  of  the  chain, 
o 


ck  ■  ‘  r 
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(kr2,4, . ,2N) 


and 


si 

"  ■  IT 


(2.7) 


(2.8) 


In  this  form  the  Hamiltonian  exhibits  an  explicit  dependence  on  the  mass  ratio 
u  which  is  a  convenient  form  for  our  discussion  of  the  weak  coupling  limit.  Ke 
now  transform  to  the  fully  diagonal  form  with 

9  *  ?  Vi  P  "  ^  Vi 

(2.9) 

q'  -  L'  xkjhj  pj  •  £’  (k»l  ,3, . .  2N-1) 

^2k  “  n2k  *  p2k  "  n2k  (k"l »2»  •  •  •  »N) 
where  the  prime  on  the  sums  indicates  j«0,l ,3, . .2N-1 . 

By  direct  substitution,  one  then  finds  that  the  secular  equation  may  be 
written  in  the  form^ 


G(z)  -  0  -  z  -  n2  -  Tn(z) 


(2.10) 


with 

V*> 


(2.11) 


2 

The  solutions  (z  -  s£)  of  this  secular  equation  are  the  eigenfrequencies  (s^ 
the  fully  diagonalized  Hamiltonian.  In  addition  one  finds  that 


*ok  *  nr*  oo 
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with 


1  ^  ■>  >  r 

n  «  |  »•  *  \  • 

i«i  (sk-u.) 


Aw  indicated  in  the  introduction  we  shall  focus 
correlation  function  p  of  the  impurity .  The  normalized 
given  by 


(2.13). 

on  the  momentum  auto¬ 
correlation  function  is 


P(0 


<rr(t)> 

<p‘> 


(2.14) 


whe-.e  the  bracket  denotes  the  canonical  average 

<A(M.pV)>  ■  (2.15). 

JdPdQdp  dq  c  P 

If  one  now  uses  the  equations  of  motion  of  the  harmonic  lattice  and  the 
indicated  canonical  transformation*  [Eqs.  (2.1)  through  (2.13)]  one  finds  that 


PM(b.t)  ■  X*  cos  s.t  (2.16) 

"  k  **  *  (k»0,l , . . . ,2N-1) 

2 

where  X^  is  given  by  Eq.  (2.12)  and  where  the  (s^)  are  the  normal  mode 
frequencies  of  the  fully  diagonalized  liamiltonian.  The  representation  (2.16) 
is  convenient  for  discussing  the  momentum  autocorrelation  function  for  the 
finite  N  case  if  one  can  sum  the  functions  and  of  Eqs.  (2.11)  and  (2.13). 
These  sums  have  been  evaluated  and  yield  for  G(z)  and  SN(ak)  ^ 


G(ak)  ■  0  ■  (1-p)  -  ycot  3^^  |c0t1,0k  '  c,ci,°k  I 

(k«l,3, . , . ,2N-1) 


and 


W  ■ 


2cos 


2  nok 


[ 


so. 

2N  ♦  cot  sinsaj,  -  cosso^ 


] 


4  (k-l,3,...,2N-l) 

where  tho  Oj,  arc  related  to  the  normalized  eigcnfrequencies  *k  by 


(2.17) 


(2.18) 
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k  -  (1,3,...,2N-1) 


(2.19). 


2  .  2  nak 

sk  "  Sin  2(2N*i) 

It  is  easy  to  see  that  as  p+1  one  just  recovers  the  equal  mass  nearest  neighbor 
spectrum 

■  5‘"2  USTTT 

To  study  the  effect  of  the  mass  ratio  n  in  the  thermodynamic  limit, 

N-*®,  it  is  more  convenient  to  represent  the  momentum  autocorrelation  function 
by 

■  m  f  d2  u-21) 

where  the  contour  is  chosen  as  a  circle  in  the  complex  z-plane  enclosing  all 

the  zeroes  of  G(z).  The  representation  (2.21)  is  equivalent  to  that  in  (2.16). 

Note  that  in  the  limit  N->«  the  zeroes  of  G(z)  become  dense  on  the  interval 
2 

(0,uQ)  and  one  may  then  contract  the  above  contour  to  an  integral  running  just 
above  and  below  the  real  axis  in  this  interval.  If  one  computes  G(z)  from 
Eqs.  (2.10)  and  (2.11)  in  the  limit  N-"»,  one  obtains 

[**  2  1/2 

p(t)  ■  I  -K  a  cosxTdx  (2.22) 

2v  J- 1  (l-2ii)x2+p 

where  we  have  introduced  the  scaled  time 

t  »  u)Qt  (2.23) . 

This  spectral  representation  holds  for  0  <  m  1;  for  u  >  1#  a  light  impurity, 
one  obtains  another  term  corresponding  to  an  isolated  frequency  which  gives 
rise  to  a  purely  periodic  component  in  the  correlation  function. 

This  completes  the  summary  of  the  equations  for  the  oscillator  chain 
with  nearest  neighbor  interactions.  In  the  next  Section  we  discuss  the  results 
of  computer  solutions  of  Eqs.  (2.16)  and  (2.22). 
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m.  I'Al-CHI  \1  luS  Ui  AUiOOKkliUTiON  lUXtTlOJi*. 

Kc  first  consider  the  aoaentua  autocorrelation  function  of  tho  Hti  N 
particle  in  the  thcraodynaaic  liait  with  the  help  of  Kq.  (2.22).  This  integral 
cannot  he  evaluated  analytically.  For  c  I  it  lias  an  expansion  in  the  fona^’*^ 


■  V>  •  fffsf  t  <•-*«>  V’> 

This  expressions  simplifies  for  certain  anas  ratios;  for  u»l,  equal  aass  particles, 
one  obtains 


PU)  •  J0(t)  (S.2) 

and  for  n  ■  j 

p(i)  -  J#(t)  ♦  J2(0  (S.S). 

In  tho  combined  msl  coupling  liait,  irdl,  t«,  inconstant,  and  theraodynaalc 
liait ,  N«,  one  obtains  froa  the  spectral  representation,  Eq.  (2.21),  the 
veil  known  result^ 


-in  P(t)  ■  e”*1 
me 


0.4). 


The  physical  basis  of  the  utak  coupling  limit  Is  the  existence  of 
processes  occuriag  an  different  tins  scales.  Processes  which  occur  slowly 
relax  to  oquilihritm  In  the  "neon  field**  of  the  fast  processes.  In  the  present 
exanple  the  equal  uses  particles  relax  an  a  tins  scale  of  l/o#  chile  the  nansntun  of 
the  heavy  particle  (M  »  a)  relaxes  on  a  tine  scale  at  soured  by  the  aass  ratio 
a  ■  u/M.  This  latter  tins  scale  Is  slew  coopered  to  the  l/u#  tine  scale  of 
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the  mss  9  particles,  la  the  Unit  when  w«0  the  tooentue  of  the  heavy  particle 
is  a  constant  of  the  notion. 

To  obtain  infomotion  about  the  fora  of  the  itoacntua  autocot  ri-1  at  ion 
function  for  non-zero  w,  u  c  1,  we  have  nuscricafcly  integrated  l.q.  (2.22). 

In  Hfurr  1  we  plot  p(t)  versus  t  for  variou*  values  of  u.  The  tire  i* 
acasurnl  in  units  of  the  Eaxiaun  frequency  of  the  lattice,  vhtch  is  typically 
of  order  ID*15  sec.  The  danped  oscillator)*  behavior  is  clearly  evident  fren 
this  figure,  but  for  u«0.1  one  already  has  behavior  suggesting  exponential  decay. 
A  log  plot  of  this  data  in  Figure  (2)  shows  that  for  all  intents  and  purposes 
one  has  readied  the  weak  coupling  iis.it  when  A  cocparison  of  the  values 

of  eap(-ui)  and  p(i)  of  Eq.  (2.22)  for  s.«0.1  in  Table  I  shows  that  they  differ 
by  only  1%  for  i  li,  i.t.  for  tines  longer  than  about  bxlO***'  seconds. 

Note  that  there  is  an  initial  transient  period  where  the  autocorrelation 
function  aunt  ho  Gaussian  for  any  On  the  other  hand  for  long  tines,  i.c. 

tines  t  grantor  than  the  exponent ial  decay  tine  1/a.  one  has  a  correction  of 
<fi<  oscillatory  fom  to  the  weak  coupling  Unit.  Those  contributions  to 
e(t)  are  however  lots  than  n.i%  of  the  Initial  value  p(0) . 

It  should  bo  noted  that  tha  uass  ratio  appears  as  the  square  root  in 
the  Haul  1 ton! an,  Eq.  (2.5),  to  that  one  hat  exponential  relaxation  of  the 
autocorrelation  function  fnr  n  value  of  0.5  of  the  ,,snallN  paranetcr  In  the 
Manilteeiea. 

Nr  the  ant  of  finite  N  we  nust  find  the  eigenfrcquencics  of  the 
secular  Iq.  (2.10)  which  are  the  solutions  of  the  transcendental  Eq.  (2.17). 
There  am  only  Ihi  nodes  which  are  pertinent  to  this  problem  the  tcro  frequency 
node  duo  to  the  translational  invariance  of  the  lattice  and  K  nodes  arising 
fron  the  syunetric  nodes  of  the  unperturbed  lattice.  The  anti-syeectric  nodes 
hove  o  node  at  the  petition  of  tho  heavy  nass  particle  and  thus  do  not  influence 
its  dyumical  behavior.  These  N  frequencies  (s^lfor  (k«l,2,...2tt-))  an 
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c4»il)  li'jiiJ  nxH'rii'4  lly  I  run  lq».  (2.17)  and  (2.19)  and  are  used  to  obtain 
the  kith  the  In  Ip  of  Eq.  (2.12)  and  (2.13).  The  Mscntui  autocorrelation 
funf  i  ion  u*  thou  calculated  by  cari)in(  out  the  sunnation  in  Cq.  (2.16).  In 
liguri  (3)  he  plot  the  calculated  values  of  pn(w,i)  vs  t  for  various  saal 1 
value*  of  N. 

As  has  hern  pointed  out  by  Kubin^  and  others  one  expects  that  for 
tines  less  than  those  required  for  a  signal  to  propogatc  around  the  lattice, 
tin  autocorrelation  function  should  closely  approx i note  the  lb*  result.  This 
is  due  to  the  fact  that  there  is,  to  a  good  approx 1 nation,  no  way  for  the  signal 
to  know  that  it  is  travelling  in  a  finite  lattice  until  it  reaches  the 
boundaries.  For  a  signal  travelling  the  lattice  with  the  speed  of  sound,  the 
tine  i(*«.|t)  required  to  propogate  around  a  lattice  of  2N  uass  points  is  t  *w  2N. 
For  tines  t  <  2N  the  results  for  the  finite  lattice  should  thus  be  nunerically 
indistinguishable  froa  thos>  of  the  infinite  lattice.  This  is  bora  out  by 
the  calculations  presented  here.  Airthernore,  for  relatively  snail  values  of 
a,  for  which  p(t)  decays  fairly  rapdiiy,  the  finite  N  results  arc  in  excellent 
agreenent  with  the  thernodynanic  Unit  results  even  when  N  is  only  of  the 
order  of  SO.  Thus  even  in  a  havnonic  systen  where  one  has  coherent  signal 
propogatien,  snail  N  results  are  quite  accurate  in  describing  the  dynauics  of 
local  pertuibations.  One  oust  of  ce  trse  restrict  such  statenents  to  systcaa 
with  potentials  whose  range  is  anal'  coopered  to  the  site  of  the  systsn  as 
otherwise  the  boundaries  would  be  "felt"  by  the  signal  at  all  tines. 

This  "recurrence"  tine  t  **  2N  bears  no  relation  to  the  Foiacard 
recurrence  tine  which  refers  to  the  recurrence  of  a  particular  phase  point. 

These  latter  tines  are  nuch  longer  than  any  tine  scale  considered  here.^ 


IV.  jmroUWKUATlOSi  HMCTlOft  KMt  IJOTia  IU1II  L'Ali  kANul  I’.TTkACTJON 

la  this  taction  we  turn  our  attention  to  the  ease  of  a  hamunlc  lattice 
with  long  range  interaction*,  he  con»i<kr  hen*  only  the  cate  of  the  infinite 
chain. 

The  Ibailtonian  can  be  written  as 

11  "  3>f  p2  *  2*  £  l‘f  #  i  £  £  MlAl Jqj  *4,1* 

where  we  tail*  »  Q.  he  irpose  periodic  huu;iJar>  on  the  chain  and 

attune  that 


*  Ai»-ji  ■  *»  t4-,) 

Ihit  latt  condition  expresses  the  physically  realistic  assumption  that  the 
interaction  between  two  natt  points  depend*  only  on  their  reparation  in  the 
lattice,  he  take  the  interaction  to  be  of  the  fora 


*» •  s n-ot*'1  «<t<i  .  »•!.* . 

(4.5). 

Note  that  at  y»0  rue  recovers  the  nearer t  neighbor  interaction  and  os  y*l  one 
obtains  an  interaction  of  tero  strength  and  infinite  extent.  The  noraalization 
factor  (1-y)  ensures  that  the  total  potential  energy  of  the  infinite  system 
rraains  constant  for  all  v. 

by  procedures  analogous  to  those  used  in  obtaining  Eq.  (2.22)  (  the 
spectral  representation  of  the  aoaentuo  autocorrelation  function  of  the  heavy 
mss  particle  is  found  to  be 

pCv.v.O  •  n  f  2/  3 - 3 - 

^-1  lw  ♦2w(l-w)ui(wtY)swi(l-|i)a(Ai(wtY)si2(wtY) 

where 


COSfc'tdw 

(d.d) 


-  12 


(4.5) 


•  - VTjl1  ^  .■  1  •  jt 

u-o,/*Ki*>r-4iw2] 


fcU.t)  * 


4>*j 

- y-i - - 


(4.6) 


and  khcrc  the  f rtt^uency  *•  is  nm  dincosionless  having  him  Haled  by  the 
Kiiinui  frequmwy  w#1  The  autocorrelation  function  p(p,Y,0  of  Lq.  (4.4) 
can  easily  be  evaluated  for  different  value*  of  u  and  y  by  miner!  cal  integration 
and  aooc  representative  result*  are  displayed  in  Figure*  (4).  The  *aae 
da  aped  oscillatory  behavior  a*  shorn  in  Fig.  (1)  for  nearest  neighbor  inter* 
actions  is  obtained  in  this  case  for  intemediate  value*  of  u  and  For 
stMll  enough  v  and  y  <  1  one  again  obtains  an  exponential  decay  of  the  auto- 
correlation  function.  As  v  increases,  l.c.  as  the  range  of  Interaction  increase!:, 
one  nust  go  to  snaller  values  of  n  to  obtain  exponential  relaxation.  In  the 
uc#ib  coupling  Unit,  ubich  one  can  obtain  fra*  Lq.  (4.4),  by  setting 


«'  ■  m/m  ,  y'  ■  mi 


and  koldiag  <'  constant,  one  finds 

^  M 

lia  •  (£!)•  ** 


MT«c,N«* 


(4.«). 


The  weak  coupling  Unit  thus  yields  again  an  exponential  relaxation  for  the 
MMcntua  autocorrelation  function.  For  y*0»  lq.  (4.1)  reduces  to  the  nearest 
neighbor  Interaction  result  (5.4). 

It  has  been  denonstratedl^10^  that  for  interactions  which  are  of  the 
fora  Ajj  ■  exponential  relaxation  is  obtained  in  the  Unit  w-*0  when 

the  squared  frequency  distribution  of  the  equal  mass  systca  (e*l)  satisfies 


-  IS  - 


(4.9) 


C(y2) 

for  Hill  «.  Our  interaction  Mtrii  utiifio  this  criteria  since  on«.  finds 
from  fcqs.  (4.4)  through  (4.6)  that  for 


C(u) 


2I  .ii^Q 


|(l*t)2-4Ts^Jll-w2J,/- 


(4.10) 
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\  .  MCNAHY  AM*  ni^OJSSlOX 

Starting  with  the  Hamiltonian  for  a  linear  harmonic  chain  of  2N 
I'jriido  with  one  tail  defect  particle,  we  have  calculated  via  analytical 
dynanic*  the  noteentuu  autocorrclat ion  function  of  the  tan  defect  particle  for 
finite  chain*  and  for  non -zero  nas*  ratio*  w.  he  liave  shown  that  one  obtain* 
to  a  very  good  approximation  an  exponential  relaxation  of  the  momentum  auto* 
correlation  function  for  mass  ratios  as  high  as  n«O.I  and  for  chains  with  only 
So  particles.  As  is  well  known,  passage  to  the  thermodynamic  and  weak  coupling 
limits  yields  tho  result  that  the  dynamics  of  the  infinite  mass  particle  is 
rigorously  described  by  a  l‘okkcr*Planck  equation.  This  in  turn  yields  the 
rigorous  result  that  the  momentum  autocorrelation  function  of  the  infinite  mass 
particle  has  an  exponential  tine  decay,  it  is  clear  from  tho  above  results  that, 
at  least  for  tho  harmonic  chain  studied  here,  the  Fokker*Pianck  equation  (or, 
oquivaiently,  tho  corresponding  Langovin  aquation)  can  be  used  to  describe  to 
a  very  good  approximation  tho  dynamics  of  a  heavy  (but  not  necessarily  infinitely 
heavy)  particle  in  a  fiaite  chain.  Tho  thermodynamic  limit  and  the  woak  coupling 
limit,  while  necessary  to  obtain  rltoraus  analytical  results  for  the  validity 

of  stochastic  equations  of  notions,  art  thus  unnecessarily  stringent  conditions 
for  tho  use  oi  stochastic  squat  loss  is  describing  tho  dynamics  of  the  model 
considered  here. 

One  important  question  which  lModluttly  arises  is  how  applicable  this 
conclusion  is  to  other  systems.  Is  this  true  in  general  or  are  those  findings 
quite  specific  to  the  harmonic  chalnf  Wo  believe  that  our  result  on  the 
validity  (in  an  approximate  rather  than  rigorous  sense)  of  stochastic  equations 
such  u  the  master  equation,  tho  Fokker-Plsnck  equation,  or  the  Langevin  equation 
far  outside  (whatever  that  may  mean  la  any  given  cate)  the  thermodynamic  and 
weak  coupling  limits  is  a  very  general  one.  Unfortunately,  this  must  remain  a 
conjecture  for  the  time  being  since  we  know  of  no  general  proof. 


•  IS  - 


A  somewhat  related  study  has  recently  been  carried  out  by  Bishop  and 
berne^)  vho  investigated  via  computer  calculations  the  on^et  of  Brownian 
■otion  in  a  one  dimensional  fluid.  They  found  an  exponential  relaxation  of 
the  velocity  autocorrelation  function  for  clusters  with  mass  M  *  25n  where  m 
is  the  aass  single  fluid  particle.  Their  results  thus  point  definitely  in  the 
sane  direction  as  ours. 

The  analytic  treatments  of  this  problem  imply  that  the  momentum  auto¬ 
correlation  filiation  can  be  written  as  a  series  expansion  in  the  coupling 
parameter  n  where  the  lending  term  is  given  by  the  thermodynamic  and  weak 
coupling  limit  result  and  where  the  correction  terms,  arising  from  finite 
strength  coupling  and  the  finite  size  of  the  system,  are  proportional  to  powers 
of  p .  A*  we  have  seen,  these  correction  terms  are  very  small  in  the  example 
considered  hero.  It  should  be  pointed  out  that  it  is  exceedingly  difficult  to 
demonstrate  analytically  that  corrections  to  the  "rigorous"  stochastic  equations 
are  indeed  proportional  to  powers  of  the  coupling  parameter.  It  is  omy  very 
recently  that  it  has  been  shown'  '  that  the  Langevin  equation  and  the  Pokker- 
Planck  equation  for  a  heavy  particle  in  a  classical  fluid  do  '  ve  correction 
terns  proportional  to  the  coupling  parameter  u  which  arc  bounded  for  all  times. 
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TABLE  1 


Comparison  of  exp(-pt)  and  p(p,;)  of  Eq.  (2.2 2)  for  p=0.1 


i _ cxp(-in) _ p(p,  0 


0 

1.000 

1.000 

3 

0.741 

0.819 

6 

0.549 

0.570 

9 

0.407 

0.415 

12 

0.301 

0.292 

15 

0.223 

0.212 

18 

0.165 

0.159 

21 

0.122 

0.108 

24 

0.090 

0.076 

27 

0.067 

0.056 

30 

0.050 

0.039 

33 

0.037 

0.028 

36 

0.027 

0.020 
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FHIURI:  CAPTIONS 


I  1  gu  re  2 


Figure  3 


Figure  4 


Normalized  momentum  autocorrelation  function  p(i,p)  for  p=.05, 

.1  and  .f».  The  scaled  time  i  is  measured  in  units  of  the  maximum 
frequency  c  of  the  lattice. 

A  logarithmic  plot  of  p(i,u)  vs  \  for  p=0.5  and  p=.l.  The  circles 
are  data  points;  this  data  is  closely  approximated  l)y  a  straight 
line . 

Normalized  momentum  autocorrelation  function  p^(t,p)  for  2N+1 
particles  with  N-12  and  N=24  and  p=.l. 

Normalized  momentum  autocorrelation  function  p(t,p,y)  for  u=.l  and 
Y® . 1  and  y=.5.  For  small  values  of  y  one  still  obtains  exponential 
relaxation  for  small  p. 
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